Quantum walks as massless Dirac Fermions in curved Space-Time 
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A particular family of time- and space-dependent discrete-time quantum walks is considered in one 
dimensional physical space. The continuous limit of these walks is defined through a new procedure 
and computed in full detail. In this limit, the walks coincide with the propagation of a massless 
Dirac fermion in an arbitrary gravitational field. Finally, the limiting procedure and the main result 
itself are carefully discussed. 
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I. INTRODUCTION 

The first quantum walk was built by Feynman [3fj| | as 
a possible discretization of the standard, massive Dirac 
dynamics in flat space-time. General discrete-time quan- 
tum walks have then been introduced in the physics liter- 
ature by [1] and [23] and the continuous-time version first 
appeared in [l7j . Quantum walks are the simplest formal 
analogues of classical random walks and are important in 
many fields, ranging from fundamental quantum physics 
flfjl to quantum algorithmics B, l22j |. solid state ph- 
syicsjl, @-(8| and biophysics [ITI.[l6|. 

Quantum walks have been realized experimentally, for 
example as transport of trapped ions [H, Hil j of photons 
in wave guide lattices |23 or optical networks [1^ and of 
atoms in optical lattices [20[ . Quantum walk experiments 
of two photons [25| | have recently been performed, with 
the possibility of simulating Bose or Fermi statistics [13] , 
and cavity QED quantum walks have also been proposed 
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Following Feynman's idea, several authors have stud- 
ied the continuous limit of general quantum walks. Most 
publications 0, H ©, GJ, EJ GHH] onl y 

envisage quan- 
tum walks with constant coefficients. The continuous 
limit of quantum walks with time-and space-dependent 
coefficients has been considered only recently, in [l~3l — 
Il5j . These references present several families of quantum 
walks, in both (1 + 1) or (1 + 2) space-time dimensions, 
whose continuous limit is described by a flat space-time 
Dirac equation with generalized mass term and electro- 
magnetic coupling. The electromagnetic field is gener- 
ated by the space-time dependence of the angles defining 
the walks and thus vanishes if these angles are constant. 

The other gauge field wich couples naturally to a Dirac 
spinor is evidently gravity. Yet, quantum walks whose 
continuous limit are described by Dirac equations in 
curved space-time have remained elusive. This Letter 
considers a particular family of discrete-time quantum 
walks with non constant angles in (1 + 1) space-time di- 
mensions and associates to this family a continuous limit 
which is described by a massless Dirac equation in curved 
space-time. This result opens the way to possible labora- 
tory experiments simulating the propagation of fermions 
in relativistic gravitational fields. It also establishes a 
connection between general relativity and all the afore 
mentioned fields, where quantum walks are useful. 



The material is organized as follows. The family of 
quantum walks considered in this Letter is presented in 
Section ITT] The retained limit procedure is presented and 
implemented in Section [III! The differential equation de- 
scribing this limit is identified as a Dirac equation in 
curved space-time in Section IIVI and the result is finally 
discussed in Section IVl 



II. FUNDAMENTALS 

We consider quantum walks defined over discrete time 
and discrete one dimensional space, driven by time- 
and space-dependent quantum coins acting on a two- 
dimensional Hilbert space H. The walks are defined 
by the following finite difference equations, valid for all 
(j, m) e N x Z: 
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The index j labels instants and the index m labels spatial 
points. For each instant j and each spatial point to, the 
wave function ^ j m — ^ m &L + V>5„&ii has two compo- 
nents tpf m and ipfm on the s P m basis (bz, &r) and these 
code for the probability amplitudes of the particle jump- 
ing towards the left or towards the right. Note that the 
spin basis is interpreted as being independent of j and 
to. The total probability = E m (I i>j, m ? + I ^> ?) 
is independent of j i.e. conserved by the walk. The set 
of angles {#j,m> (j, m) e N x Z} defines the walks and is 
arbitrary. 

Consider now, for all (n,j) £ N 2 , the collection Wj 1 = 
{^nj,m)me_z- This collection represents the state of the 
quantum walk at 'time' nj. For any given n, the collec- 
tion 5™ = (W")jgN thus represents the entire history of 
the quantum walk observed through a stroboscope of 'pe- 
riod' n. The evolution equations for S n are those linking 
W™ + i to Wp for all j. These can be deduced from the 
original evolution equations ([1]) of the walk, which also 
coincide with the evolution equations of S . In particu- 
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lar, the evolution equations of S 2 read: 

' l Pj+2,m — c j+l,m ( c j,m+li^j^ n +2 + s '] '.m+lV^m) 

+s j+ i, m (sj, m -i^j; m - Cj, m -i^jf m _ 2 ) , (3) 

^j+2,m — s J+l,m ( c J,m+lV^m+2 + s j,m+l' t Pf,m) 

-Cj+l, m (Sj, m -Xlpf >m ~ Cj l7n -i^ m _ 2 ) , (4) 

where Cj m = cos(0j TO ) and s jm = sin(0j m ). 



where the subscript T (resp. X) indicates a derivative 
with respect to T = t/r (resp. X = x/X). 

In terms of the spinor \1/ = ip L bL + ip R bR, these equa- 
tions read simply 



* f + (cos9)P^ x = QW, 



(8) 



where the operators M and R are represented, in the 
base (6l,6r), by the matrices: 



III. CONTINUOUS LIMIT 



P 



— cos f z sin ( 
— isin# cos 6 



(9) 



To investigate the continuous limit of S n , we first in- 
troduce a time step At and a space step Ax. We then 
introduce, for F — and F — 8, a, function F defined on 
M + x M. such that the number Fj m is the value taken by 



and 



F at the space-time point (tj — jAt, x n 
example, equation ([1} then reads: 



mAx). For 



At, x r , 

At, X r 



= B{9{t 3 ,x m )) 



i> L (tj,x m + Ax) 
ip R (tj,x m - Ax) 



(5) 

where the tilde has been dropped on functions to sim- 
plify the notation. We now suppose, that the functions 
can be chosen at least C 2 in both space and time vari- 
ables for all sufficiently small values of At and Ax. The 
formal continuous limit of S n is defined as the couple 
of differential equations obtained from the discrete-time 
evolution equations defining S n by letting both At and 
Ax tend to zero. Let us therefore introduce a time-scale 
r, a length-scale A, an infinitesimal e and write At = re 
and Air = Xe. The continuous limit of S n can then be in- 
vestigated by Taylor expanding in powers of e the discrete 
equations defining S n . For the limit to exist, all zeroth 
order terms must identically cancel each other and the 
differential equation describing the imit is then obtained 
by equating to zero the non identically vanishing, lowest 
order contribution. Consider first S , which is identical 
to the original walk. It is rather obvious that zeroth or- 
der terms cancel each other in ([5]) only if the operator 
B defining the walk tends to unity as e tends to (see 
[] for a detailled discussion of this point). The operator 
B defining the family of walks considered in this Letter 
does not depend on e and is different from unity for all 
values of 9. Thus, S 1 does not admit a continuous limit 
for the family of walks defined by |T]). But S 2 on the 
other hand does. Indeed a straightforward computation 
delivers the following equations obeyed in this limit by 
-0 L and -0 fl : 

V>t - (cos 2 6»)V4 + ^Oin 2 ^^! = 
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IV. DIRAC EQUATION IN CURVED 
SPACE-TIME 

The operator P is self-adjoint and its eigenvalues are 
— 1 and +1. Two eigenvectors associated to these eigen- 
values are 



b- = i I cos — ) 



b + =i[ sin - ) b L 



^n- | b R , 



cos - ) b R . 
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The family (£>_,6 + ) forms an orthonormal basis of the 
two dimensional spin Hilbert space. Let us now rewrite 
equation ((5J) in components, but in this new orthonormal 
basis. A tedious but straightforward computation leads 
to: 

Ox 
2 
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(cos0)i/j x 
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where ip~ and if>~ are the components of \I/ in the new 
basis. This form of the equations makes it easy to check 
that the continuous dynamics conserves the total proba- 
bility tt(T) = f x dX | * |2= J x dx (| r |2 + | ^ | 2)j as 

it should. 

Suppose now, to make the discussion definite, that 
cos 9 is non negative and introduce in space-time 
{(T,X)} the Lorentzian metric g defined by its covari- 
ant componenents 



1 







— cos t 



(14) 

defines the 

V=gdX = 

V cos 9 dX in physical ID A"-space, where g is the de- 
terminant of the metric components Dirac spinors 



where (/i, v) € 
canonical, scalar 



{T,X} 2 



This metric 



: volume' element T) g X 
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are normalized to unity with respect to T> g X, whereas 
is normalized to unity with respect to dX. We thus 
introduce $ = ^(cos^) -1 / 4 and rewrite the equations of 
motion p3[) in terms of We obtain: 



7 



z v a 



0, (15) 



where fi E {T, X}, a € {0,1}. The usual 2D gamma 
matrices are: 



7° = 



1 

1 



1 

-1 



(16) 



(17) 



and the are the components of the diad (orthonormal 
basis) eo = &t and e\ = Vcos9 ex on the original coor- 
dinate basis (ex, ex)- Equation (fT5|) is the standard (3l| 
equation of motion for a massless Dirac spinor in (1 + 1) 
dimensional space-time with metric g. 

Given any coordinate system, the metric describing 
a arbitrary gravitational field in a (1 + 1) dimensional 
space-time can be represented by the matrix of its compo- 
nents g M „ in the so-called coordinate basis. Since this ma- 
trix is symmetric, there are only three independent com- 
ponents. But these three components do not represent 
three degress of freedom for the gravitational field. In- 
deed, relativistic gravitation is invariant under a change 
of gauge i.e. a change of space-time coordinates; thus 
space-time coordinates can be chosen in an arbitrary way. 
Since there are two space-time coordinates, the choice of 
gauge can be encoded in two of the three independent 
metric components, and the remaining metric compo- 
nent then represents the single, true degree of freedom 
of relativistic gravitation in (1 + 1) dimensional space- 
times. In particular, in regions of space-time where there 
is a non-null time-like field, the metric components of 
any (1 + 1) dimensional gravitational field can be put 
under the form (fT4")l by a suitable coordinate choice. In 
other words, the single angle 9{t, x) is enough to describe 
any (1 + 1) gravitational field in regions where 'there is 
a time'. This excludes for example the horizons of black 
holes, but includes the domains outside and inside the 
horizons. 



V. DISCUSSION 

We have considered a particular family of quantum 
walks in one-dimensional physical space and introduced 



a new procedure which gives a continuous limit to walks 
from this family. We have computed this limit and proven 
that it coincides with the propagation of a massless Dirac 
fcrmion in an arbitrary gravitational field. 

Previous work on the continuous limit of quantum 
walks does not investigate the possibility of not keeping 
every step of a walk to build its limit, yet this seems to be 
the key ingredient to generate an apparent gravitational 
field. Indeed, retaining all steps for quantum walks on 
the line necessarily leads to Dirac equations in flat space- 
time, even for walks defined by time-and space-dependent 
angles [] ; the effect of the time- and space dependance is 
then, not to generate a gravitational field, but to gener- 
ate an electric field. Also, keeping all steps to define the 
continuous limit makes it impossible to define this limit 
for all walks [l3l[l4j. For example, the family considered 
in this Letter only has a continuous limit if one retains 
one step out of every n — 2k, k <E N (the case investi- 
gated here corresponds to n = 2, k = 1) . As a bonus to 
then being able to define a limit, one obtains a new and 
interesting limit i.e. a massless Dirac fermion coupled to 
a gravitational field. 

The work presented in this Letter should naturally be 
extended in several directions. One should first investi- 
gate systematically all quantum walks on the line defined 
with a quantum coin acting on a 2D Hilbert space, but 
try and build limits by keeping only one step out of every 
n > 2. One should also extend the main result of this 
Letter to quantum walks defined on physical space of 
higher dimension and/or defined by quantum coins act- 
ing on a higher dimensional Hilbert space. A particular 
goal would be to obtain walks whose limits are described 
by a Dirac fermion coupled to both a gravitational and 
an electromagnetic field. 

The geometry appearing in the continuous limit pre- 
sented in this Letter is an intrinsic property of the walk 
itself, not of the 'real', 'physical' space-time which might 
be used to realize the walk experimentally. Another ex- 
tension of this work would therefore be to consider walks 
defined on curved physical spaces (resp. graphs), and 
to investigate how the geometry of the underlying space 
(resp. graph) couples to the intrinsic geometry of the 
walk. This is not a purely academical problem, since 
quantum walks can model photon transport in networks 
of algae, which may have a non trivial geometry. 

The main result of this Letter also suggests that con- 
cepts from general relativity and diferential geometry 
may play a key role to understanding the behaviour of 
quantum walks on graphs and their use in quatum al- 
gorithmics. This role should also be investigated thor- 
oughly. 
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